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ABSTRACT 

A set of linearized equations of motion to predict the 
linearized dynamic response of a single rotor helicopter in a 
hover trim condition to cyclic pitch control inputs is described. 
The equations of motion assume four fuselage degrees of freedom; 
lateral and longitudinal translation, roll angle, pitch angle, 
four rotor degrees of freedom; flapping (lateral and longitudinal 
tilt of the tip path plane), lagging (lateral and longitudinal 
displacement of the rotor plane center of mass) and dynamic 
inflow (harmonic components). These ten degrees of freedom 
correspond to a system with eighteen dynamic states. In addition 
to examination of the full system dynamics, the computer code 
supplied with this report permits the examination of various 
reduced order models. The code is presented in a specific form 
such that the dynamic response of a helicopter in flight can be 
investigated. As described in the report, wdth minor 
modifications to the code the dynamics of a rotor mounted on a 
flexible support can also be studied. 


INTRODUCTION 


The equations of motion for the linearized small 
perturbation motion of a single rotor helicopter about a hovering 
trim condition were formulated using a Lagrangian approach. For 
a linearized investigation of the hovering flight dynamics it can 
be assumed that the fuselage center of mass remains in an 
horizontal plane during the disturbed motion if collective inputs 
are not of interest. Vertical translation will be uncoupled from 
horizontal translation and fuselage roll and pitch. Thus the 
collective pitch is assumed constant at its trim value. Also as 
a consequence of this assumption, the coning angle and the 
average value of the induced velocity are also assumed to be 
constant in the code. The influence of the tail rotor is 
neglected and the yaw angle is assumed to be constant at its trim 
value, as coupling of the degree of freedom will be relatively 
weak near hover. Also the steady value of the lag angle is 
assumed constant at its trim value. The documentation of a 
larger code for hover and translational flight trim conditions 
which includes the vertical and yaw degrees of freedom as well as 
the tail rotor effects is nearly complete and will be available 
in the near future. 

The fuselage center of mass is assumed to lie on the rotor 
shaft and no fuselage aerodynamics are included. Thus the hover 
trim condition corresponds to zero values of longitudinal and 
lateral cyclic pitch. Only the specific formulation for a 
helicopter in flight is considered in this discussion. See the 
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description of modifications to the code 
e other problems such as ground resonance, 
rotor dynamics are included in the code through the use 
lade coordinates. The longitudinal and lateral tilt of 
plane relative to the shaft are described by the 
es a^, b lg . The lag motion is described by and 

ding to lateral and longitudinal displacement of the 
ter of mass. These are the first four state variable in 
ter program. Four body coordinates are fuselage pitch 
elage roll (d> r ), lateral translation (y ) and 

r t 

nal translation (x^). Note that longitudinal 

on is positive to the rear. This particular choice of 

e depends on the selection of the elements in the T 

described below. Two state variables v and v 

c s 

the dynamic inflow. Thus the displacement variables for 
system are 

{a ls' b ls' »V r 2- V V y F ’ X F ’ V V 


The dynamic inflow variables are first order while the other 
variables are second order so that there are a total of eighteen 
state variables. The system equations are ordered such that the 
rotor and body displacement coordinates are followed by the rotor 
and body velocities and then the dynamic inflow components. Thus 
for the full system the state variables appear in the following 
order in the computer code output, 
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V 


V 


{*i> T = { a ls » b is> y i« y 2 • e F’ *F’ y F ’ X F ’ a ls’ b ls’ r l* r 2’ 

0 F ’ ^F ’ y F ’ X F ’ V c’ V s } 

The control variables are, 

<V T = < A ls' B ls> 

For the reduced order system option without dynamic inflow, 

the last two state variables (v ,v ) are dropped and there are 

c s 

sixteen state variables. 

For the quasi static option, both the rotor blade motion and 
the dynamic inflow become algebraic variables and are eliminated 
from the equations of motion to yield eight state variables, 
describing the fuselage motion 

- • • " T 

{ X q S } = { 0 jr> ’ y F’ X F’ ®F' ^F ’ y F’ X F^ 

The quasi-static case has two options, with and without dynamic 
inflow. 

Note that since there is no dependence in the equations of 
motion on fuselage displacement (y^, x p) so * n these cases 

there will be two zero eigenvalues. 

While the quasi-static model is useful for developing 

physical insight, it has been shown in Reference 1 by comparison 

of this theory with flight test that the rotor-body coupling and 

dynamic inflow are significant in predicting the actual response 

of a helicopter to cyclic pitch. 

The equation of motion are formulated for an articulated 

rotor helicopter with equal flap and lag hinge offset. Hingeless 
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rotor helicopter dynamics can be approximated by the addition of 
springs about the flap and lag hinges. These spring constants 
appear in the input file. 

The rotor blade element aerodynamics are assumed linear, 
i.e., the blade element lift coefficient is assumed to be a 
linear function of blade element angle of attack and the blade 
profile drag coefficient is assume constant. With these 
assumptions the blade aerodynamic forces can be integrated along 
the blade span to the flapping hinge and multiblade coordinates 
introduced to describe the rotor motion. The dynamic inflow 
model used in the computer code is described in Appendix I and is 
due to Peters. Rotor angular velocity is assumed constant. 
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NOMENCLATURE 


Rotor Geometry 

0 n = 9 - A, siniji - B, cosiJj + A9„ + A0,, 

B o Is Is H E 

ill = rotor blade azimuth angle measured from longitudinal axis 
(x F ) 

0 q = collective pitch assumed constant, calculated in code 
( TNOT ) , rad. 

A ls ,B ls = lateral and longitudinal cyclic pitch, control input 

terms corresponding to lateral and longitudinal stick 
deflection, radians 

A8^ = blade pitch change due to hinge coupling 

A8 h = d* Ap + e* A<f 

d* corresponds to a 6 ^ hinge 
e* corresponds to hinge 


A0 e 

= blade pitch 
plate 

change due to 

elastic deformation 

A0 e 

= A0 cos4* + 
c 

A0 sincjj 
s 



> 

CD 

O 

= (A - 1) 0x H 

+ B0y H + 

Cx H + 

Dy H 

A0 

s 

= E0x„ + (F - 
ri 

l) ey H + 

g* h + 

Hy H 


For rigid 

shaft A = 

1, F 

= 1, B=C=D=E=G=H=0 


P = a q cosip - b^ s sintK flap angle, positive up. rad 


coning angle assumed constant, calculated in code (BETA), 
rad 


a, , b, = multiblade flapping coordinates (state variables) a, 
Isis 1 » 

longitudinal tilt, positive for flap back, rad 
b i ^ lateral tilt, positive for tilt down on right 
side, rad 

i = i - y, cos - Yo sinij), lag angle, positive for lag opposite 
o 1 c 

to rotation, rad 

i = steady lag angle assumed constant, calculated in code 
o 

(RNOT), rad. This quantity does not appear elsewhere in 
the code due to the use of multiblade coordinates 

^1 ’ r 2 ” ^ultiblade lag coordinates corresponds to lateral 

displacement of the rotor center of mass, positive to the 

right, and correspond to longitudinal displacement, 

positive forward, rad 

Inflow 

v = v + v x cosiJj + v x simjj, positive downard, fps 
o c s 

v q = average induced velocity, assumed constant, calculated in 
code ( VNNT ) , fps 

v , v = harmonic components, tip amplitude, calculated as 
c s 

described in Appendix I 

Fuselage 

h = distance between fuselage CG and hub, ft 
Hub Motion 

x„ = longitudinal displacement of hub, positive aft, ft 

y = lateral displacement of hub, positive to the right, ft 
H 

0Yu ” pitch angle, positive nose up, rad 
h 
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0x = roll angle, positive left side down, rad 
H 

Generalized Coordinates 


X H 

= z 

tx i q i 

y H 

= z 

ty i q i 

0y H 

= z 

X q i 

e*H 

= z 

‘ex. q i 

i 

t X . 
1 

, ty 

i’ t 0y i ’ 


j. 

generalized coordinates (q^) to hub 
motion. These are selected in the code as 

q = 0 , fuselage pitch angle, positive nose up, rad 

1 F 

q = 4) , fuselage roll angle, positive right side down, rad 

2 F 

q = y , fuselage center of mass lateral displacement, 

3 F 

positive to the right, ft. 

q = x , fuselage center of mass longitudinal displacement, 

4 F 

positive to the rear, ft 

T = {tx., ty., t 0x , t Qy }, transformation matrix 


h 

0 

0 

1 


0 

h 

-1 

0 


0 

1 

0 

0 


Note that the sign of the roll angle coordinate has been reversed 


in the transformation matrix. 



DISCUSSION 


In order to investigate various coupled rotor-body problems 

such as: 1) The dynamics of a helicopter disturbed from hovering 

flight; 2) The dynamics of a helicopter resting on the ground 

(ground resonance); 3) The dynamics of a rotor on a flexible 

mount, the effect of hub motion is treated in the following way. 

Hub motion is assumed to be horizontal displacement in two 

directions (x,y), and rotation in pitch and roll ( 9 , 9 ), 

y x 

(x H ) = (x, y, e y , 0 X ) T - 

Generalized coordinates (q^) are used to describe the support 

motion and are related to the hub motion by a transformation 
matrix T: 

(x H > = [T] {q.} 

The matrix T can be selected to represent a variety of supports. 
The equations of motion are of the general form: 

Rotor Equations: 

a rr* x r } + a rh {x h^ = 0 

Support Equations: 

a srr {x r^ + a shr {x h* + a shs* x h* = 0 

The A 1 s are second order operators with the following 
interpretation 


RR = r0t0r 

forces 

and 

moments 

due 

to 

rotor motion 

= rotor 

forces 

and 

moments 

due 

to 

hub motion 




A = support forces and moments from rotor due to rotor 
SRR 

mot ion 

support forces and moments from support due to hub 
motion 

support forces and moments from rotor due to hub 
mot ion 


SHS 


SHR 


Now to express the equations of motion in terms of support motion 


{q^} the transformation 


{x H } = [T] {q.} 

is introduced and the equations of motion expressed as: 




+ < A RH T){q i l 1 ° 


^SRR^R 1 


+ (i SHS T + A SHR T1 * q i* " ° 


This is the general form of the coupled rotor-support equations 
of motion. Thus, the equations of motion of a rotor on a rigid 
support with no hub motion (q^ = 0) are: 

a RR {x R> = 0 

and the equations of motion of the support without the rotor are 

A SHS T {q i^ = ° 


It is assumed in the equations of motion given below that the 
matrix of operators describing the support 

system motion is diagonal and is of the form: 

a shs t = \ [S ii s2 + 6 ii s + R li’ 

Thus the q.’s are generalized coordinates describing the support 
i 

degrees of freedom* Care must be taken in the selection of the 
support motion coordinates, such that the ma Crix (A shs T) is 
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diagonal. The factor 1/2 appears due to the normalization of 
* 

of support equations of motion. 


1 

2 


M. . 


1 1 


bl, 


(M. . ) 

li 


The system matrices are of the general form, 


M 


RR 


M 


[M] 


M 


SRR 


RH 


(M SHS + m shr ) 


etc . 

Potentia l Ener gy and Virtual Work Due to Thrust 

The potential energy terms and the virtual work due to the 
rotor thrust must be treated carefully since second order terms 
are involved which are sensitive to the assumptions made 
regarding support motion. The rotor terms are calculated based 
on the assumption that the hub motion takes place in horizontal 
plane. The linearized equations of motion assuming either 
horizontal motion of the fuselage support or horizontal motion 
the hub, are the same as shown below. The potential energy and 
virtual work due to thrust must be treated carefully. These 
terms will be different in the two cases but their sum will be 
the same. It is necessary to include in the evaluation of thes 
terms the effect of a vertical displacement (z). It is desirab 
to consider the formulation with the fuselage or support center 
of mass translational displacement in a horizontal plane. This 
will also be a suitable approximation if the rotor is on a 


flexible mount. 


Rigid Fuselage/Shaft - Rotorcraft in Flight 


First consider the case in which the equations of motion are 
used to examine the dynamics of a helicopter in flight where the 
shaft is assumed to be rigid. Both vertical motion of the hub 
and fuselage are permitted at this point. The potential energy 
of the system is (Fig. 1): 

V = M f g z F + bMgg z H (1) 

and the virtual work due to the thrust is: 

= T cos9 sine dx u - T sine dy 
T x y H x H 

( 2 ) 

+ T cosG cosG dz„ 
x y H 

The thrust, T, is assumed constant since the trim condition is 
hovering flight, however since it has a non-zero value in trim 
the virtual work due to the second order displacement dz^ must be 

included since it will produce a first order term in the virtual 
work . 

The relationship between the fuselage displacements and the 
hub displacements are: 


X H = 

X F 

+ 

h’ 

cosG sin0 
x y 


V- 

SC 

n 

y F 

- 

h’ 

s in0 

x 

(3) 

N 

SC 

II 

Z F 

+ 

h’ 

COS0 COS0 

X y 


height 

of 

the 

rotor 

center of mass 

above the fuselage 


center of mas s . 

S B 

h’ = (h + B ) 

m b ° 

It can be shown that the linearized equations of motion that 



result from assuming either hub motion constrained to a 
horizontal plane (z^ = 0) or fuselage motion constrained to a 
horizontal plane (z p = 0) will be the same if the potential 

energy and virtual work terms are evaluated consistent with 

either of these approximations. The vertical velocity of the hub 

(z„) does not appear in the rotor kinetic energy terms due to the 
H. 

hub motion assumption used in deriving the rotor terms in the 
equations of motion. Using the alternate assumption that the 
fuselage center of mass motion is constrained to the horizontal 
plane (z = 0), there is no linearized contribution to the 

r 

kinetic energy from z which arises only due to rotation of the 

H 

shaft. Therefore, it is only in the potential energy and virtual 
work terms which arise due to the thrust vector that the effect 
of this different assumption must be considered. The effect of 
rotation on vertical displacement must be included in the 
evaluation of these terms. 


Horizontal Fuselage Motion 

First consider the evaluation of the potential energy and 
virtual work based on the assumption that the fuselage center of 


mass remains in a horizontal plane (Zp =0). In this case, 


potential energy terms will arise due to the vertical 

displacement of the rotor mass. 

Equation (1) becomes, 

V = bM n g z„ = bM gh’cosG cos0 
b H d x y 


14 


ae ae 

— = ( -bM„ gh ’ ) ( s in0 cos0 - - + cos8 sin© • " ) 

aq i B x y aq t x y aq i 


s - bM D gh ’ (0 — — — + 0 — — *-) 
B x aq t y dq i / 


These potential energy terms are considered apparent spring terms 
in the computer code since they depend upon the support model 
assumption. They are entered directly in the input file as 
terms. They take the following form, 

1^— = - bM gh ’ ( ( £ t. .q.) t. . + (I t Q .q.) t. .) 
aq i B 0X1^1 0X1 0yi l 0yi 

For the coordinate selection given below, 

aq x = 00^ = ” bM B gh ’ 0 F 

aV - dV - vm ■ 

- a*; - ’ bM B gh *f 

The virtual work terms due to rotor thrust simplify for this 
assumption of (z r = 0) by noting that since the rotor thrust lies 

along the shaft, angular motion of the shaft does not contribute 
to the virtual work. Using Equations (3), with = 0, equation 

r 

(2) becomes, 


6W_ = T cos 8 sin0 6x„ - T sin0 6y„ 

T x y F x J F 


Q 6x„ + Q 6 y r 

x F F y F 


Consequently the linearized terms are 
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equat ion 


x 


F 


= t e 


= T 


I t n . q. 
0yi i 


Yp equation 


(5) 


T = * T 1 ‘«xi "i 


These terms are entered directly in the spring matrix [K]^ for 
for the specific form of the T matrix given below, i.e., 


1 0y 1 1 


t 0x2 1 

In the computer code these are terms RK1(8,5) and RK1(7,6). 

The transformation matrix is selected in this case to 
identify the as follows: 


q x = < 0 F ’ V y F ’ X F }T 

where 

0^ fuselage pitch, positive nose up 

r 

<t> fuselage roll, positive right side down 
F 

y fuselage lateral translation, positive to the 
right 

x r fuselage longitudinal translation, positive to the 
F 

rear 


and the matrix T is 




T. . 
i J 


T = 


[t 


xi 


h 

0 

0 

1 . 0 


yi 

o 

h 

1.0 

0 


0xi 

0 

- 1.0 

0 

0 


t 0yi ^ 


1.0 

0 

0 

0 


Note that the sign convention for the roll angle has been 
reversed (t^ 3 = “ 1.0). 

With this selection of coordinates the support terms are of 
the form 

( M . . q . + K . . q . ) 

11 1 111 

specifically 


q i <v 


dy S F - Kg V 


q 2 (* F ) 


*F - K B V 


q 3 <y F> 


< M F v f> 

q 4 ( V 

X F^ 

where Kg is obtained from the potential energy terms (equations 
( 4 ) ) above . 


K b = - bMg gh’ 
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The spring constant K D is negative, representing the fact that 

D 

the rotor mass is above the fuselage mass. 

Horizontal Hub Motion 

v It will now be shown that if the alternate assumption 

(z„ = 0) is employed the same terms will result. In this case 
H 

equations (1) and (2) become: 

v V-MpgZp 

6W m = T cos© sin© 6x tI - T sin0 6 y„ 

T x y H x H 


and from Eq. (3) 


v 

6z F 

= + h’ 

s in 0 

cos 9 
X y 

60 

X 

+ h 1 c o s 0 sin© 60 

x y y 



6x h 

= 6Xp - 

h’ 

sin8 sin© 
x y 

60 + h’ cos 0 c o s 0 

x x y 

>> 

CD 

t O 

V 

6y H 

= 6y F ' 

h’ 

COS0 60 

X 

X 




Assuming small 

angles 

and 

comb in ing 




6V - 

6W t = 

(Mpg 

- T) h’ 

0 

X 

60 

X 


w 


+ 


- T) h’ 

0 

y 

60 

y 

(6) 



+ 

T 0 

y 

6x F - T 

e 

X 

6y F 



From vertical equilibrium 

T = M p g + b Mgg (7) 


it can be seen that equations (6) and (7) give the same result 
for the two terms given by equations (4) and (5). The terms are 
6V - 6 W t = -bMggh’ 9 x 60 x 

- bMg gh* 0 y 60 y (8) 

+ T 0 y 6Xp - T 0 x 6y F 
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Rotor Mounted on Support - Fixed Support Point 

Now consider the case in which the rotor is mounted on a 
flexible transmission/shaft system. There is assumed to be no 
vertical motion of the support point (z^ = 0). The potential 
energy is due only to the blade mass, 

V = bMgg Z H 


The virtual work due to the thrust is given by equation (2). If 
we assume that the flexibility of the mounting is described by a 
number of normal modes of amplitude q^, then the hub motion is 
related to the q^ by the transformation matrix T. The support 
equation of motion would then be of the form, 


M. . {qV } + C . . {q . } + K . . { q . } + {f^ Q. } 

11 i li i li i dq . i 


+ ^RR^R 1 


* a shr t { V = 0 


where the effect of hub vertical motion (z u ) should be included 

rl 


in the evaluation of the potential energy term due to the blade 

mass and the virtual work term due to the thrust. In general as 

noted above it is not necessary in the linearized formulation to 
include the effect of z^ in the remaining terms due to the rotor. 

^ii’ ^ii’ ^ii are 2 e neral ized mass, damping and stiffness 

associated with the ith mode. In general z„ will be a nonlinear 

H 


function of shaft deflection, i.e., a function of the lateral 
deflection of the shaft squared. Physically this means that the 
z^ deflection terms associated with thrust component and the 


potential energy terra will give rise to apparent spring terms. 



Rigid Shaft 




V 


Consider first the rigid shaft case as discussed above. For 


small deflections, 

6W t = T 9 y 6 x h - T 9 y 6y H + T 

6V = bMgh ’ 6 z r 


( 9 ) 


For rigid rotation of the shaft with no translation of the 
support motion ( dx^ = 0, = 0) 

6x h = h’ 60 y 

6y H = -h ’ 60 x CIO) 

6z H = - h ’ 9 x 69 x - h ' 9 y 69 y 


Substitution of equations (10) into (9) shows that the virtual 
work is zero. An apparent spring is produced due to the 
potential energy term associated with blade mass. That is, 
physically, if the shaft is rigid the thrust does not produce a 
moment about the support. Thus to examine this problem the terms 
corresponding to the virtual work in the computer code should be 
removed (RK1(8,5), RK1(7,6)). 

Flexible Shaft 

Consider now a simple case in which the shaft has one 
deflection mode, in longitudinal plane only. The shaft 
deflection is represented by a mode shape t}^. The local 
deflection of the shaft is (Figure 2): 

x s = ^(Z) q x 

The local slope is 
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v' 






d-n, 

e = = XT - ( z ) 
s dz I 

At the hub location, the mode shapes are represented by the T 
matrix coefficients. 

Til (h) = T 


'1 

dn- 


xl 


-(h) 


dz * 9y 1 

The z deflection of the hub is 

L 


6z 


H 


where 


ax 

az 


Ii 7<5 s ~ ( & 2 dz > 

1 o 


d-n . 


Thus 


dz 


L dn 


1 % 2 


SZ H = - < f <d^> dZ > -i iq i 
o 

This is an apparent spring terms or an effective change in 
stiffness due to the potential energy and virtual work terms. 

Note that the thrust is a "follower" force so that the 
apparent change in stiffness depends upon the mode shape. The 
virtual work due to the thrust is ( 6y^ = 0) 

T y H H 


Ther ef ore , 


«* T * T [T eyl - / (g^-) 2 dz] q, 6q, 


diK L di)| 9 

* T ldF- (L) V L) - s <dT-> d2 > q l dq 

o 


'1 


1 
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For a rigid shaft 


z 

11 1 = L 

and 

6W t = ° 

For a flexible shaft the sign of the apparent spring depends upon 
the sign of the term in brackets above. This term can be also 
expressed as, 

T A 2 

L d 7). 

<5W t £ T ( J* n 1 g dz ] q i 6q i 

o dz 

^ by integrating by parts in the case Tij(O) = 0* 

Thus the thrust will often tend to act as a negative spring 

even though a tension is produced due to the fact that it is a 

" "follower" force (Figure 2). 

These terms will be small, for a typical flexible shaft on 

which a rotor is mounted. That is, the thrust and blade weight 

~ will not have a strong effect on the shaft frequencies. 

In the general case, with more complex shaft /transmiss ion 

modes, the T coefficients correspond to the deflection and slope 

w at the end of the shaft or support. If the mode characteristics 

are only available numerically then an approximation must be made 

to determine the z u dependence with deflection which should be of 

H 

the general form indicated above. Note also that if the shaft is 

assumed rigid with a rotation about the base, the above equations 

for 6 z„ will recover the correct result for a rigid shaft, 
ti 
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The generalized stiffness and mass are: 


K n = f EI(n ") 2 dz 
o 


M 11 = ^ dz 
o 
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SUMMARY 


A linearized set of equations of motion for a coupled rotor- 
body system for the trim condition of hovering flight have been 
described. The rotor kinetic energy terms are derived based on 
the assumption the hub motion takes place in a horizontal plane. 
The formulation permits the analysis of various coupled rotor- 
body problems, including a rotor mounted on various supports and 
free flight. Support motion can be assumed to lie in a 
horizontal plane such that vertical displacement of the hub 
arises only from rotation (rigid shaft) or support flexibility if 
the potential energy and virtual work due to the thrust force 
account for the vertical displacement of the shaft. For this 
assumption, vertical displacement of the hub will be of second 
order and will not contribute to the kinetic energy, but will 
give rise to first order term in the potential energy and 
virtual work due to thrust. Other virtual work terms due to 
inplane forces and hub moments are not affected by this 
assumpt ion . 

The terms to be evaluated for the specific support 
assumption are: 

V = M f k z f + bMgg z H (1) 

and 

6W t = T cos© x sinQy 6x R - T sinG^ <5y H 

( 2 ) 

+ T cos9 cos0 6z„ 

x y H 

The equations of motion as described in detail below and 
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currently in the computer program assume rigid shaft motion and 

fuselage motion in a horizontal plane (z r = 0). This assumption 

r 


yields the relationships in equation (8) 


6V = bM n gh ’ (- 9 60 - 0 69 ) 

B x x y y 


( 8 ) 


6W t = T 0 y 6 x f - T0 x 6y p 

The terms due to the virtual work of the thrust are entered 
directly in the matrices (RK1(8,5), RK1(7,6)) and their placement 
is reflected by selection of the q ^ ’ s as noted. The potential 

energy terms give rise to apparent spring terms that are entered 
in the input file Again their placement is reflected 

by the selection of the q^’s. 

Thus to study dynamic problems of a rotor on a support with 
other selections of generalized coordinates, the potential energy 
and virtual work due to the thrust must be calculated for the 
specific problem and entered suitably in the computer program. 

The virtual work terms are in a form corresponding to free flight 
and are due to thrust. In general they should be removed from the 
computer code (RK1(8,5), RK1(7,6)) # for other problems. 
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PROGRAM INPUT FILE AND OPTIONS 


The computer program is set up to provide the equations of 
motion for a hovering helicopter in four cases: 

1. Full system with dynamic inflow 

2. Full system without dynamic inflow 

3. Quasi-static system with dynamic inflow 

4. Quasi-static system without dynamic inflow 

The full system includes blade flap and lag degrees of freedom, 
body degrees of freedom and dynamic inflow. The dynamic inflow 
model is described in Appendix I of this report and the quasi- 
static formulation is given in Appendix II. For the full system 
with dynamic inflow, the state variables are, 

{x} = {a x , b r r r r 2 , q L , q 2 , q 3 , q 4 , v c , v s } T 

The selection of the T matrix yields for the q’s, 


q l 0 F 


q 2 ~ 


q 3 = y F 


q 4 = X F 


For the full system without dynamic inflow, v c and v g are dropped 
as state variables. For the quasi-static system, 

{X qs } = {q l’ q 2 ’ q 3 * q 4 }T 


The rotor blade motion and the dynamic inflow are not state 
variables and are eliminated as explained in Appendix II. Note 
however that in the quasi-static case there is the option to 
include quasi-static inflow terms which will of course influence 


the results. 

The data input file describes the physical characteristics 
of the helicopter. The sample included is for the case of a 
helicopter in flight with a rigid shaft. 

The first three lines of the data file are the generalized 
mass, spring and damping associated with the q. degrees of 

freedom. These depend upon the choice of the T matrix (the next 
two lines of the data file). For the selection of the T matrix 
elements in the sample data file, 


q l 

II 

CD 

'TJ 

, M n 

= I 

yy 

q 2 

= * F 

’ M 22 

= i 

XX 

q 3 

= y F 

’ M 33 

= ra F 

q 4 

X F 

• M 44 

nip 


Note that these are inertial characteristics of the helicopter 
without rotor blades. The generalized spring terms arise from 
considerations discussed above and are equal to 


K 11 K 22 " bM B g (h + Mg *o> 


The generalized damping terms are equal to zero in this case. 
The next two lines are the elements in the T matrix, 


t x 


ty. 


1 0x . , 1 0y . 

1 1 


The selection in the example is explained in the introduction. 
The sixth line in the data file contains 


V 
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mechanical lag damping f t-lb/rad/sec 


= equivalent lag spring ft-lb/rad 

= equivalent flap spring ft-lb/rad 

= rotor rpm, rad/sec 
The seventh line is, 

M_ = blade mass, slugs 

D 

M = blade first mass moment, slug-ft 
s 

2 

I = blade flapping moment of inertia, slug-ft 
B 

(flap and lag moments of inertia assumed equal) 

The eighth line is, 

R = blade radius, ft 
e = flap/lag hinge offset, ft 
c = blade chord, ft 

a = blade solidity (number of blades appears through 
solidity) 

a = blade lift curve slope 

The ninth line is, 

3 

p = ambient air density, slug/ft 
6 = blade average profile drag coefficient 

The tenth line is, 

d* = pitch change due to flap (6^ hinge) 

e* = pitch change due to lag 

The eleventh and twelfth lines permit introduction of swash plate 
deflections with shaft bending. In the form given there is no 


K 


HL 


K 


HB 


n 


swash plate deflection relative to the shaft with shaft motion. 
The thirteenth line is: 


T = trim thrust, lbs 

h = dynamic inflow parameter 

WRF = wake rigidity factor 

h is the height of the equivalent cylinder of air to be 
accelerated and WRF is the wake rigidity factor which should be 
set equal to 2.0 as explained in Appendix I.. 
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APPENDIX I 


DYNAMIC INFLOW 

The equations of motion also permit addition of dynamic 
inflow components. Since the trim condition is hovering flight 
only variable harmonic components are introduced. The harmonic 
components are assumed to vary linearly with radius. 

Av = v xcosiJj + v x simp 
c s 


The equations describing the harmonic inflow component 
amplitudes at the blade tip are: 

4C m 

7 t V + V = - k [ ] 

I c c I ao 


4c a 

t t v + v = - k T [ — ] 

I s s I ao 


The harmonic components are positive for downward flow 

through the rotor. ti is the time constant associated with 

response time of the inflow to changes in aerodynamic moments on 
the rotor (C^, C^) and kj is the proportionality between the 

aerodynamic moments and the inflow in the steady state. These 
are expressed as 

h 

T I = 

2v [) f 
o w 

. aoRft 

k I = 

2v f 
o w 

There are two variable parameters in these expressions whic 
are input parameters to the computer program 


_ 

h = p where H is the height of an equivalent cylinder 

of air which must be accelerated. The theoretical 
value is H = .453R. 


w 


is the wake rigidity factor. The correct value of 

this parameter is 2.0 (non-rigid wake) although in 

some places in the literature it is taken as 1.0 

(rigid-wake). This terminology is from Miller. 

(f = ( WRF ) ) . 
w 

In addition the harmonic induced velocity terms will change 

the aerodynamic forces on the rotor blades. Denoting v = {v^, 

, T 

V s } • 

The complete equations of motion with dynamic inflow are: 


M 

0 


0 

0 



(x| K DYE ~1 

V v { + l -DYB 1 -DYC ] [ v ) 


B 

DF 

Introducing the notation 

4C. 


MB 


MB 


ao 




where C 


MB 


the nondimens ional blade flapping moment has two parts, 

.* 


one due to the rotor and body motion, denoted C 
to the harmonic inflow. Therefore, 

.* 


MBo 


and one due 


dC 


'MB 


= C 


MB - 


MBo 


dv 


Note that in the equations above v is dimensional 
notation of the computer program. 

.* 


In the 


DYE v 


rn 


2 dC 


MB - 
v 


dv 


rfi 

2R 


ac 


MB 


dv 
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where 




DYC (1,1) = ( VNNT ) ( WRF ) 2 + £j> (RMFVC + EB FNVC) 

HH i (K 


DYB = - PP C UD = DYB2 x + DYB 1 x 

MB 

O 


PP = 


aaRfi^ ^1 


WRF 


w 


VNNT 


average induced velocity, dimensional 


HH = h = h R, dimensional 


h and WRF are input constants in computer program 


ac 


MB 


av 

ac 


c 

* 

MH 


- - RMFVC 


= - [RMFVC + ( EB ) ( FNVC ) ] 


av 


For the quasi-static inflow assumption, 

DYCv = - DYB 

and v can be eliminated from the force and moment equations, 
additional terms in the force and moment equations accounting 
inflow effects are, 

- ( DYE ) (DYC ) _1 DYB 
The term due to v^ equals, 


The 

for 


v 
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(-PP) 


rfi r aC MB 

2R av } 


1_ A L ( ^MH 

t t t t flH ' - 
II 3v 


) 


-} c 


MB 


Substituting for PP 


rfi 


^1 aC MB 

nR - 

av 


1 , k I 8C MH 

rr - 


dv 


-} c 


MB 


nR 


ao 


2 v f 
o w 


ac 


MB 


av 


ac 


MH 


av 


- RMFVC = |- 


rn 


a a 

8v f 
o w 


1 + { 


a a 


8v f 
o w 


“} c 


MB 


This term can be combined with other aerodynamic terms in the 
equations of motion which are of the form 

rft 2 c * 

2 MB 


to yield, 


r {- 


l + 


ao 


8v f 
o w 


2 

-} 2- C* 

} 2 U MB 


Thus the effect of the inflow in the quasi-static case can 
be viewed as a reduction in Lock number. The effective Lock 


V-/ 


33 


V--' 


number is , 


w 


r E 


1 + 


ao 

8v f 
o w 
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APPENDIX II 

QUASI-STATIC FORMULATION 

The quasi-static formulation is obtained by first separating 
the rotor and body degrees of freedom. 


{a x , b 1? Yj, Y 2 ) 


V = {v , V } 
k c s 


X b = {q r q 2 ’ q 3 ' q 4 } 


The equations of motion are written as: 


Mn M 12 


M 21 M 22 



C 11 C 1 2 


C 21 C 22 



K 11 K 12 


K 21 K 22 



+ [DYE] {v} = 


{u} 


The equations of motion for the inflow components are: 

{v} = [ DYC ] { v } + [ D YB 1 ] { x } + {DYB2}{x} + [DF] {u} 

The quasi-static assumption is obtained by setting: 


M 


11 


M 21 ‘ C 11 ' C 21 ' 0 


v = 0 


solving the inflow equation, 

{v} = - [DYC] -1 [ D YB 1 ] { x } - [DYC] -1 [DYB2]{x} - [DYC] -1 [DF] {u} 
This is substituted into the first set of equations: 

[M] {x} + [C-DYE DYC -1 DYB2]{x} + [K-DYE DYC -1 DYB 1 ] {x} 


[F + DYE DYC 1 DF] {u} 
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Rewriting the equations: 


m {x} + [C] {x} + [Rj {x} = [P] {u} 


or 


ft ll fi 12 


^21 ^22 


C 11 e i2 


e 21 6 22 



Ki K 12 


K R 

i K 21 *22 



r V 


L F 2 J 


{u} 


The quasi-static assumption is: 

«11 - "21 ■ e il * e 21 ’ 0 

Solving the first equation for the rotor motion, 

{ x r } = - Kjj fl 12 (* b ) “ Kjj <x b ) - R bl Rj 2 {x b ) 

+ RjJ ?! {u} 

This is substituted into the second equation, 

f ^22 “ ^ 21 ^11 ^ 12 -* ^ X b ^ + ^22 ~ ®"21 ^11 ^ 12 ^ ^ X b * 

+ [R 22 - K 21 r i;i r 12 i {x b } 


= f f 2 ' K 21 K il < u > 


These are the equations of motion for the quasi-static system 
[QM] {x } + [QC] {x } + [QK] {X } = [OF] {u} 

These matrices are calculated in the subroutine QUASI. 


COMPUTER PROGRAM 


Subroutine Matrix 2 calculates the following matrices 
[AX] [BX] [CX] [FX] 

[DF] [DYB] [DYC] [DYE] 

The complete system equations are: 

[AX] {x} + [B] (x) + [CX] (x) + [DYE] {v} = {FX} u 

{v} = [DYC] {v} + + [DYB 1 ] {x} + [DYB2 ] {x} + [DF] {u} 

THE VARIABLES ARE 

T 

{ ^ } — { & ^ i > ^ 2 5 ^ i * ^2 * q 3 ’ ^ 

{V} = (v c , v s } 

NOTE: 

AX, BX, CX, FX are renamed M, C, K, F 
when Matrix 2 is called. 

Complete system eq ua tions of motion 


* 


0 1 0 

n 


' 0 1 

• 1 

X 


-M _1 K -M _1 C -M _1 DYE 

i 

. 

X 

) + 

M _1 F 

* , 


DYB 1 DYB2 DYC 


l V 


. DF 
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SAMPLE DATA FILE 
(BHEFA.DAT) 


UH-60A BLACKHAWK PARAMETERS 




38512.0, 4659.0, 460.9, 460.9 

-7959. 0, -7959. 0, 0. 0, 0. 0 

0 . 0 , 0 . 0 , 0 . 0 , 0.0 

6.87, 0.0, 0.0, 1.0, 0.0, 6.87, 1.0, 0.0 

0 . 0 , - 1 . 0 , 0 . 0 , 0 . 0 , 1 . 0 , '? 0 < 0-0 
4600, 0. 0, 0. 0, 27 . 0 
7.98, 86, .70, 1512.6 

£6.83, 1.25, 1.73, 0.0821, 5.73 


1. 

95E 

—03, 

0. 

01 

5 

0. 

0, 

0. 0 




1 . 

0, 

0.0, 

0. 

0 , 

0. 0 

0. 

0, 

1. 0, 

0. 

0 

, 0. 0 

15370 

. 0, 0 

. 

* , £ 

. 00 


BLRCKHfiWK data for flight test. Includes two additional parameter 
following the thrusts HH, the height of the equivalent cylinder o 
air; and UlRF, the wake rigidity factor, equal to 1.0 for a rigid 
ar id elm 0 for nori—r i g i d wake. 


.f 

w a k e 





V 


EIGENVALUE 

I 


4 

nr 

s 

7 

a 

9 

^ 10 

1 1 
IE 

13 

14 

15 

W 1 tj 

1 7 
18 


□PEN-LDDP EIGENVALUES 3- 


REAL 

I MAG 

. 000E+00 

. 000E+00 

. 000E+00 

. 000E+00 

~9. 095E+00 

5. 203E+0 1 

-9. 095E+00 

-5. 203E+01 

-1 . 9Q3E+00 

3. 91 1E+01 

-1. 983E+00 

-3. 91 1E+01 

— 57GE+01 

£. 464E+00 

— £’. 576E + 01 

-£. 4G4E+00 

— 1 . 353E + 00 

1 . 82QE+01 

-1 . 353E+00 

-1. 8E8E+01 

-£. 997E+00 

4. 940E+00 

— £'. 997E+00 

-4. 940E+00 

— 4. £'63E + 00 

. 000E+00 

-1. 51 1E+00 

. 000E+00 

5. 173E-02 

3. S75E-01 

5. 173E-02 

-3. 275E-01 

6. 505E-03 

3. 539E-01 

6. 505E-03 

-3. 539E-01 




SAMPLE PROGRAM OUTPUT 
FULL SYSTEM WITH DYNAMIC INFLOW 
INPUT DATA FILE (BHEFA.DAT) 

f 

E igenva 

lue Identification 


(1,2) 

Associated with lack of dependence on x f , y ? 


(3,4) 

Advancing flap mode 

V 

(5,6) 

Advancing lag mode 


(7,8) 

Inflow mode 


(9,10) 

Regressing lag mode 


(11,12) 

Coupled roll/body flap mode 


(13,14) 

Coupled pitch/body flap modes 


(15,16, 

17,18) Coupled roll, pitch, body translation modes 
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Notes on the usage and compilation of H5NEW: 

H5NEW is an updated version of a program for calculation of helicopter rotor/body linearized dynamic 
equations of motion. The helicopter is assumed to be in hover; formulation is sufficiently general to 
accommodate a wide variety of helicopter and rotorlypcs, as explained in the detailed documentation. 
These paragraphs are meant to aid installation of the program on a machine other than that on which it 
was developed, namely, an IBM-PC with an 8087 numeric coprocessor. 

Files on this distribution diskette include source code for the generation of matrices in the equations of 
motion, as well as executable code that should run, unmodified, on an IBM-PC class machine with an 
8087 math chip installed. Scenarios are presented below if this is not your current installation. 

If you will be running the program on another machine, the source must be recompiled and relinked into 
an executable module. One of the notable problems that will be immediately encountered is that the 
source for the eigenvalue/eigenvector routines is NOT on this diskette, due to space constraints. These 
routines, called from the subroutine "EICSYS", are part of the standard EISPACK matrix subroutine 
package for real general matrices, and should (hopefully) be available at your installation. Note that the 
eigenanalysis portion of the main program H5NEW is really a post-processing function, and not intrinsic 
to the generation of the system matrices in the first place. Thus, one could just comment out the lines in 
the main program to eliminate the eigenanalysis, and then recompile and relink the code. 

The file "LNK" contains a "batch"-like file for the linking process, following the successful recompilation 
of the program and subroutines on the disk. Linking of the object code ("'.OBJ) can be done via. 

link @LNK 

provided that the EISPACK routines have been compiled and grouped into their own library. 

Compilation of the programs on another computer (under a UNIX or VMS environment) will of course 
depend upon the appropriate incantations necessary for the resident FORTRAN compiler, and are best 
formulated by a quick study of the contents of the listings. 

The main routine is called "H5NEW", which reads in both controlling information and names of input 
data files. This main executive then passes control to "MATRIX2”, which is given the task of computing 
all of the system matrices for the problem. Then, according to the wishes of the user, one or more support 
routines are used to combine the matrices produced from ''MATRIX2". These include: FRST, a 
subroutine that reduces the mass-spring-damping formulation into 2-N sets of first-order ordinary 
differential equations; MINV, a matrix inversion routine, required in several matrix combination 
operations; QUASI, a subroutine that generates a "quasi-static” model based upon the full rotor+body 
dynamics equations; and EIGSYS, a subroutine that in turn calls EISPACK routines to compute system 
eigenvalues and eigenvectors for either the full or quasi-static systems. 

Also on the disk is a file called OPTSYS.EXE, which is an executable program (again, requiring an 8087 
chip) for computing system poles for matrix-vector representations of linear systems. This routine allows 
one to investigate effects of certain feedback on the pole locations of the closed-loop system. 

BHEFA.DAT is an example data file for input to the program, with some additional commenting after the 
last line of numeric input. 

Should installation questions arise, you may contact me for further help: 


Bob McKillip 

Mechanical and Aerospace Engineering Dept. 

Princeton University 

P.O. Box CN5263 

Princeton, NJ 08544-5263 

(609) 258-5147 


